Spectra of Open-Shell Nuclei with Pade-Resummed Degenerate Perturbation Theory 
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We apply degenerate many-body perturbation theory at high orders for the ab-initio description of ground 
states and excitation spectra of open-shell nuclei using soft realistic nucleon-nucleon interactions. We derive a 
recursive formulation of standard degenerate many-body perturbation theory that enables us to evaluate order- 
by-order perturbative energy and state corrections up to the 30th order. We study 6,7 Li as test cases using a 
similarity renormalization group (SRG) evolved nucleon-nucleon interaction from chiral effective field theory. 
The simple perturbation series exhibits a strong, often oscillatory divergence, as was observed previously for 
ground states of closed-shell nuclei. Even for very soft interactions resulting from SRG evolutions up to large 
flow parameter, i.e. low momentum scales, the perturbation series still diverges. However, a resummation of the 
perturbation series via Pade approximants yields very stable and converged ground and excited-state energies in 
very good agreement with exact no-core shell-model calculations for the same model space. 

PACS numbers: 21.60.De, 21.10.-k, 02.70.-c 



I. INTRODUCTION 

The precise description of the nuclear spectroscopy is one 
of the major challenges in nuclear structure theory. A vari- 
ety of many-body approaches, both, exact and approximate, 
have been developed to tackle this challenge. A simple and 
well-known tool to obtain approximate solutions of the many- 
body problem is Rayleigh-Schrodinger many-body perturba- 
tion theory (MBPT) Qj. Low-order MBPT has been used for 
studies of systematics of ground-state properties of closed- 
shell nuclei throughout the nuclear chart IJ-Q], as well as for 
infinite neutron and nuclear matter calculations |Hl @]. Re- 
cently, we have used high-order MBPT, i.e. an order-by-order 
evaluation of the perturbation series up to typically 30th order, 
to address ground-state energies of light closed-shell nuclei 
0]. It turns out that the perturbation series in general does not 
converge. However, through a resummation by Pade approx- 
imants we have shown that one can utilize the information of 
the divergent power series to determine the ground-state en- 
ergy of closed-shell nuclei with the same precision as in the 
no-core shell-model (NCSM) approach B8I 410I1 using the same 
model space. 

Unfortunately, standard nondegenerate Rayleigh- 
Schrodinger perturbation theory is applicable only to 
ground states of closed-shell nuclei. Both, the step to excited 
states and the step to open-shell nuclei leads to degeneracies 
in the unperturbed (zeroth-order) energy level, which must be 
considered explicitly. In this paper, we investigate degenerate 
Rayleigh-Schrodinger perturbation theory to overcome this 
limitation. We derive recursive formulas for an efficient 
order-by-order construction of the perturbation series. With 
this extension we are able to study ground and excited states 
of closed and open-shell nuclei. We focus on the application 
of degenerate many-body perturbation theory (DMBPT) up 
to high orders and Pade resummations for the description 
of ground states and excitation spectra of light nuclei, 



specifically 6 Li and 7 Li, and we compare the DMBPT results 
with exact NCSM calculations for the same model space 
and Hamiltonian. We employ a nucleon-nucleon interaction 
from chiral effective field theory fill [l2tl at next-to-next-to- 
next-to-leading-order (N 3 LO) in the version of Entem and 
Machleidt 11311 after an additional similarity renormalization 
group (SRG) transformation [Til [l5ll to further soften the 
interaction. As expected the DMBPT series does not 

converge, not even for very soft interactions. However, we 
will show that a resummation through Pade approximants 
leads to stable and accurate predictions for the spectra. 

This paper is organized as follows: In Sec. [TTJ we present 
the formalism of DMBPT, derive the relevant formulas for 
the energy and state corrections, and highlight their recursive 
structure. The latter allows us to construct the perturbation 
series up to high orders. In Sec.|III]we introduce the Pade ap- 
proximants that we use to resum the typically divergent power 
series from DMBPT. Finally, we show results for the spec- 
tra of 6 Li and 7 Li obtained from the DMBPT power series in 
Sec. 1171 



II. DEGENERATE MANY-BODY PERTURBATION 
THEORY 

The starting point for the formulation of MBPT is the eigen- 
value problem of the intrinsic Hamiltonian 



H\V n ) = (T-T om + V) PP„> = £JY n >, 



(1) 



with r cm denoting the center-of-mass kinetic energy and T - 
T cm being the intrinsic kinetic energy. Note that the interac- 
tion V is general and might include three-body forces as well. 
Next, we partition the Hamiltonian H into an unperturbed part 
Ho and a perturbation W with an auxiliary parameter A such 
that the original Hamiltonian is recovered for A — 1 



H 



partitioning 



H A =H + AW 

= Ho+A(H-H ). 



(2) 
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Formally one has complete freedom in defining this partition- 
ing and, thus, in choosing the unperturbed Hamiltonian Hq. 



2 



In practical applications the choice of Hq is often motivated 
by computational simplicity. The eigenvalue problem of the 
unperturbed Hamiltonian 



H |<&„> = e n |<D„) 



(3) 



<0W|, yielding 

g„-i 



J] ( (®nd' I W |O nc > - E^6 d , e ) X nd,ne = , 



(8) 



e=0 



defines the unperturbed basis { |<t>„)} that is used through- 
out the perturbative expansion. This basis should be suf- 
ficiently easy to handle formally and computationally, but 
also adequate for the physical system under consideration. 
For nuclear-structure applications typical choices for Ho 
are Hartree-Fock or harmonic-oscillator (HO) single-particle 
Hamiltonians. Throughout this paper we choose the latter, 
i.e., the unperturbed states { |<1>„)} are given by Slater determi- 
nants of single-particle HO states. The energy eigenvalues e„ 
are determined by the sum of the single-particle energies of 
occupied states. 

In the simplest form of MBPT the unperturbed state that 
represents the eigenstate of interest is required to be nonde- 
generate. When using the HO basis, this is true only for the 
ground states of light doubly-magic nuclei, e.g. 4 He, ie O and 
40 Ca 0]. However, for open-shell nuclei or excited states of 
closed-shell nuclei the unperturbed states exhibit degenera- 
cies and one has to resort to degenerate Rayleigh-Schrodinger 
many -body perturbation theory. In the following, we derive an 
iterative formulation of DMBPT that enables us to efficiently 
evaluate perturbative corrections to the energies and states up 
to very high orders. 

To characterize the degeneracy of the unperturbed states, 
we introduce an additional degeneracy index d for the un- 
perturbed Slater determinants |<t> m /), which labels the states 
spanning the g„ -dimensional degenerate subspace associated 
with the energy eigenvalue e„. As a consequence, the index 
d also appears in the power-series for the perturbed energies 
and states 



nd nd 



E„d(A) - Si 



(4) 
(5) 



with E^j = e„. We insert this ansatz in the eigenvalue prob- 
lem of the Hamiltonian (0 and match same orders of A. For 
order A we recover the unperturbed eigenvalue problem (0, 
whereas for order A p with p > 1 we obtain 



where we used the orthogonality (dWlCw) = &dd'- Equa- 
tion ([8]) is an eigenvalue equation in the degenerate subspace 
for e„. The eigenvectors define the expansion coefficients 
Xnd.ne of the unperturbed states in Eq. ([7), and the eigenvalues 
the first-order energy corrections E\]. Moreover, the follow- 
ing relations hold for the unperturbed states | V F (I ^) 



Xnd,ne Xnd'.ne ~ &dd' 

e=0 



(9) 



Using the intermediate normalization ( > P^| V F„ £ /) = 1 we ob- 
tain from Eq. (0) the relation 



(10) 



W'X) = So P ■ 

After multiplying Eq. (0) with ( V P^ | and using the previous 
orthogonality relation, we obtain a simple expression for the 
p-th order energy correction 



7 (p) 

"nd 



(11) 



which has the same form as in nondegenerate MBPT J3]. The 
derivation of the perturbative corrections to the states is 
more involved. We start off with formally expanding the p-th 
order state corrections for p > 1 in terms of the unperturbed 
basis. Within the degenerate subspace « of the target state we 
use the unperturbed basis \ y f!j^) of Eq. (0, for the orthogo- 
nal subspaces m + n we use the naive unperturbed basis of 
Slater determinants \<f> me ) for simplicity. Thus we obtain the 
following expansion of the p-th order perturbative correction 



rr&n 



(12) 

The remaining task is to derive expressions for the expansion 
coefficients 



W|^- 1) )+Hol^ ) ) = S^I ,I t" J; >- 



u(P-j)\ 



7=0 



(6) 



The unperturbed states |*P?5) obviously enter in Eq. (0). How- 
ever, for each g„ -dimensional degenerate subspace for an un- 
perturbed energy e„ we can choose arbitrary linear combina- 
tions of the naive Slater determinants |O m /) to represent the 
unperturbed states 



s„-i 



1*0 = J]Xnd,ne\®ne), 
e=0 



(7) 



In order to fix the expansion coefficients Xnd,ne we consider 
Eq. (0 for p = 1, insert the expansion (0) and multiply with 



C 



(p) 

ndjne 



vOmrf) m + n, 



and 



D 



(p) 

nd,ne 



= <¥«M> e + d. 



(13) 



(14) 



We multiply Eq. (0 with (3> me | for m + n, use Eq. (0 and 
obtain 



c 



ip) 

nd,me 



i 



To get the coefficients D ] we multiply Eq. (0 with (T^ 
for e + d and make use of the expansion ( fT2l and our results 



ip-j) 

nd ^ nd.me 



(15) 
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from the diagonalization in the degenerate subspace, yielding 



D 



(p) 

nd,ne 



j(i) 

^ nd 



7 m 



/ j net nd,ne 
./=' 



(16) 



Note that Eqs. (Q3]) and (dU) hold for p > 1. 

For implementation we manipulate Eqs. ( fill . (1151 1. and (fTol l 
into a more convenient form, using matrix elements of W with 
respect to the naive Slater determinants \(S?„d) as input. We 
obtain for the pth-order perturbative energy correction 



and 



g„-i g„-i 

e=0 e '=0 



(17) 



(18) 



The expressions for the C-coefficients are recast to C (p , ® = 0, 



C 



(P=i) 

nd,me 



(19) 



and 



c 



(P>2) 
nd,me 



-i) 

m' e' 



e'±d g„-\ 

- Z Z Xne',ne"(^me\ W \®ne»)D^ 
e' e"=0 



■ Y E (i ]c (p : j) 

/ i nd nd,me 



(20) 



lower-order D-coefficients and on same-order C-coefficients. 
To proceed to high orders it is indispensable to make use of 
this recursive structure, i.e. to construct the perturbation se- 
ries order by order. 

We start the construction of the perturbation series for the 
energy with the zeroth-order contribution simply given by 
the HO energy of PP^). We obtain the first correction 
from the diagonalization in the degenerate subspace. To go 
to second-order energy correction we need the first-order co- 
efficients CrJ only. To compute the energy correction of or- 
der p with p > 3 we first calculate the coefficients D^ d for 
which we need the already known C^T and all previous D- 
coefficients. Then, we use the D (p r 2) coefficients and all C- 

nd,ne 

coefficients up to order (p - 2) to calculate the coefficients 
Cf~ . According to Eq. (fT8l , the coefficients C~T enable 



'ml 



the computation of the pth-order energy correction E ( t p J. This 
scheme allows for the iterative setup of high-order DMBPT. 
We present our results based on DMBPT corrections up to 
30th order for excited states of 6 Li and 7 Li and the according 
spectra in section [TV] 

From a computational point of view, the critical operations 
that have to be performed for evaluating the perturbative cor- 
rections are matrix-vector multiplications of the Hamilton ma- 
trix and the coefficient vectors, see e.g. Eqs. < fT8l > and d20b . 
The coefficient vectors of all previous orders need to be stored. 
This is similar to a simple Lanczos algorithm lfl7ll for de- 
termining the few lowest eigenvalues and the corresponding 
eigenvectors of the many-body Hamiltonian matrix, as used 
in the NCSM. Therefore, the limitations in terms of particle 
number and model-space size for high-order DMBPT calcu- 
lations are similar to those of the NCSM. In order to apply 
DMBPT at high orders beyond the domain of NCSM-type cal- 
culations, one has to devise alternative ways of efficiently car- 
rying out the nested sums involved in the expressions above. 



III. PADE APPROXIMANTS 



For the D-coefficients we get D 
1 



(P=0) 
nd,ne 



Oand 



D 



0>i) 

nd,ne 



^ nd 



^m,e' e"=0 



nd,me' 



-J]e u : 1) d {p - 



nd.r 



(21) 



For the construction of the perturbation series up to high 
orders, carrying out all summations explicitly is very ineffi- 
cient. Closer inspection of Eqs. ( fTTb - (fJJJ) reveals their re- 
cursive structure: For the /5th-order energy correction ( fTTb we 
need the state correction of (p - l)th-order. However, we see 
from Eqs. (JT7J and ( fT8l that the energy correction depends 
only the coefficients C 7 which implicitly require the coef- 



ficients D ip d ^. The expansion coefficients C , 
depend on C-coefficients, D-coefficients and energy correc- 
tions of lower order, while the coefficients D (P } 1) depend on 



<P>D 



themselves 



Even though the recursive formulation of DMBPT will al- 
low us to perform calculations up to very high order of the 
perturbation series, simple partial summations of this series 
do not guarantee convergence of the energies. On the con- 
trary, for ground-state energies of closed-shell nuclei we have 
shown j7j] that even soft interactions lead to strongly diverging 
power series at the physical point A = 1. In this sense MBPT 
does not converge in many practically relevant cases. 

However, the coefficients of the power series of MBPT con- 
tain all the relevant physics information on the problem. An 
efficient and elegant way to extract this information are Pade 
approximants lfl8LIT9h . Instead of expanding the function E(A) 
in a simple power series, we can employ a slightly more in- 
volved ansatz and represent E(A) as a rational function of two 
power series in A 



E(A) = 



«o + Aai + A a,i . 



(22) 



bo + Ab\ + A 2 bi 
If we truncate the power series in the numerator at order L and 



4 



in the denominator at order M this defines the so-called Pade 
approximant 



[L/M](A) 



ciq + Aai + A"a2 ■ ■ ■ + A l cil 
b + Abi + A 2 b 2 ... + A M b M 



(23) 



In order to determine the unknown coefficients a, and bj from 
the known coefficients E (p} of the simple MBPT power series 
(|4|i, we use a Taylor expansion of ( 1231 and extract equations 
connecting the two sets of coefficients by matching the dif- 
ferent orders in A. Starting from a perturbation series up to 
order p this allows us to extract the coefficients of Pade ap- 
proximants with L + M < p by solving a set of coupled linear 
equations. An equivalent but more elegant way to compute the 
Pade approximant uses the ratio of determinants constructed 
directly from the corrections E ip) of the MBPT series lfl8ll 



[L/M](A) : 



£(L-M+1) 
£(L-M+2) 



£ (L) 

L-M 

2 E {p) A M+p 

p=0 



£(L-M+2) 



£(L+1) 

L-M+l 

2 E ip) A M+p - 1 

p=0 



2 E {I »A" 

p=Q 



£(L-M+1) £(L-M+2) 
g(L-M+2) g(L-M+3) 



E (L) 

A M 



£<L+l) 
A M-l 



E (L + 2) 



£(L+M) 
1 



(24) 

with E 1 - 1 '' 1 = for p < 0. Evaluating the above determinants 
for the physical point A — 1 directly provides us with a Pade 
resummed [L/M] approximation for the ground- and excited- 
state energies using the perturbative corrections E (p) up to or- 
der p - L + M. 

As we have shown in Ref. 0], Pade approximants provide 
a reliable tool to obtain the ground-state energy of doubly- 
magic nuclei in excellent agreement with exact NCSM cal- 
culations even if the simple MBPT power series exhibits a 
strongly diverging behavior. Using the recursive formulation 
of DMBPT derived above, we will extend these studies in the 
following section to ground states and excitation spectra of 
open-shell systems using the example of 6 Li and 7 Li. 



IV. SPECTRA OF 6 Li AND Li 

As a first application and benchmark of the recursive for- 
mulation of DMBPT for open-shell nuclei and excitation 
spectra, we consider 6 Li and 7 Li as test cases. In these sys- 
tems exact NCSM calculations are easily possible and serve 
as a reference to study the accuracy of the DMBPT results. 
Therefore, we perform the DMBPT calculations using the 
same many-body basis and model-space truncation as used in 
the NCSM. The unperturbed Hamiltonian Ho consists of the 
kinetic energy and a one-body harmonic oscillator potential 
leading to an unperturbed basis |<D /K /) of Slater determinants 
consisting of harmonic-oscillator single-particle states. As in 
the NCSM, we truncate the Hilbert space to a finite model 



space by imposing a maximum excitation energy Af max /zQ 
above the lowest unperturbed energy JsT IIoll . In this unper- 
turbed basis, degeneracy emerges if we study ground states of 
nuclei away from the harmonic-oscillator shell closures or if 
we investigate excited states. 

For the full Hamiltonian H given in Eq. ([TJ we include a 
two-nucleon interaction V based on chiral effective field the- 
ory at next-to-next-to-next-to leading order (N 3 LO). Starting 



from the N LO interaction of Entem and Machleidt [ 13 ] we 

" US 



apply the Similarity Renormalization Group (SRG) |Hl 
to soften the interaction by a continuous unitary transforma- 
tion. As a result, the full Hamiltonian has very favorable con- 
vergence properties as we increase the size of the many-body 
model space by increasing 7V max . Furthermore, the unitary 
transformation is believed to facilitate the order-by-order con- 
vergence of MBPT [20], which we will come back to later 
on. 

First, we investigate 6 Li in the framework of DMBPT up 
to 30th order. The degenerate subspace for the lowest unper- 
turbed energy (« = 0), which is simply the 7V max = subspace 
in the language of the NCSM, consists of 10 Slater determi- 
nants |<X>od). In a first step, we diagonalize the perturbation 
W in this subspace yielding the first-order energy corrections 



Ex] and the unperturbed basis l 1 ^)- We assign the degen- 
eracy index d — 0, ...,9 in ascending order of the first-order 
energy. Then we use the recursive formulation of DMBPT 
derived in Sec. [II] to compute the perturbative corrections to 
the energies and states up to 30th order. 

The results of the DMBPT calculations for the 6 Li ener- 
gies of all states with n — are shown in Fig. Q] for a model 
space with 7V max = 8 and fiQ. = 20MeV. We use two different 
Hamiltonians including SRG-evolved chiral NN interactions 
with SRG flow-parameters a = 0.04 fm 4 (A » 2.24 for 1 ) and 
0.16 fm 4 (A 1.58 fm -1 ), respectively. For comparison we 
show the exact NCSM results for those Hamiltonians in the 
same model space as horizontal lines. In all cases the pertur- 
bation series diverges. However, we can distinguish different 
characteristics of the partial sums of the perturbation series as 
function of the truncation order p. 

Let us first consider the harder interaction with a = 
0.04 fm 4 (blue discs). One class of states (d = 1,5,6,7,9) 
exhibits an apparent alternating convergence with increasing 
order p up to p * 12, but then at p « 16 the size of the 
perturbative corrections explodes and the partial sum diverges 
in an oscillatory pattern. For another class {d = 0, 2, 8) the 
oscillatory behavior sets in earlier and the amplitude first in- 
creases slowly before the rapid divergence sets in. A third 
class (d = 3,4) diverges monotonously starting already at 

The general situation is the same for the second Hamil- 
tonian using a SRG-transformed chiral NN interaction with 
a = 0.16 fm 4 (red diamonds). This interaction is generally 
considered to be very soft and shows a rapid convergence 
of the energies as function of model space size. Sometimes 
these interactions are termed 'perturbative', based on an anal- 
ysis of Weinberg eigenvalues in two-body systems I20I l2lll . 
Our order-by-order calculation in DMBPT up to p = 30 
shows that the softness of the interaction does not guarantee 



(0h 
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FIG. 1: (color online) Energy from DMBPT power series truncated at order p for the energy levels of 6 Li corresponding to the degenerate HO 
n = subspace using a Af max = 8 model space for Kl = 20 MeV. The indices d are determined by diagonalization in the degenerate subspace, 
i.e. by the first-order energy corrections. The dashed lines correspond to the NCSM result in the same model space. Two SRG flow parameters 
are shown: 0.04 fm 4 (A « 2.24 fnT 1 ; •, ) and the softer 0. 16 fm 4 (A « 1.58 far 1 ;♦, ). 



convergence — not even a systematic improvement of the con- 
vergence behavior. For all states we observe a strong oscil- 
latory divergence for high orders of DMBPT. Thus, in terms 
of the order-by-order perturbation theory for a light nucleus, 
even these very soft interactions lead to a divergent perturba- 
tion series and are non-perturbative in this sense. 

The divergence of the perturbation series makes it impossi- 
ble to determine a robust and unambiguous approximation for 
the exact eigenvalues from the high-order results. This prob- 
lem was already found and addressed in our previous study 
focusing on ground-state energies of closed-shell nuclei 01 • 

Even low-order estimates obtained from the second- or 
third-order calculations do not provide a reliable guideline. 
For the harder interaction with a = 0.04 fm 4 the exact 
NCSM eigenvalue is typically between the second-order and 
the third-order estimate, i.e. the second-order approximation 
gives an energy below and the third-order approximation an 
energy above the exact energy eigenvalue. However, with in- 
creasing degeneracy index d the second-order approximation 
moves up and appears above the exact eigenvalue for d — 9. 
This trend is even more pronounced for the soft interaction 
with a = 0.16 fm 4 , starting from d — 5 the second-order ap- 
proximation and all other low-order approximations are above 
the exact energy. In particular, going from second- to third- 
order DMBPT enlarges the discrepancy to the exact NCSM 
result. Thus, in general, low-order approximations do not pro- 
vide a controlled estimate for the exact eigenvalue. 

To overcome the convergence problems of the simple per- 
turbation series defined in Eq. (0), we compute the Pade ap- 
proximants according to Eq. ( f24l > and evaluate them at A = 1 . 
The diagonal [L/L], the super-diagonal [L/L+ 1] and the sub- 
diagonal [L/L - 1] Pade approximants are shown in Fig. [2] 
again for the ten states corresponding to the lowest degenerate 
subspace n — of 6 Li for the interaction with SRG parame- 
ter 0.04 fm 4 . Recall that the information that enters the [L/M] 
Pade approximant is exactly the same as in the perturbation 
series truncated at (L + M)th order. 



If we use only low-order DMBPT energy corrections as in- 
put the quality of the Pade approximants is comparable to the 
simple perturbation series truncated at these orders. If we in- 
clude information from higher orders of DMBPT, the agree- 
ment of the Pade approximants with the exact NCSM results 
improves successively. Beyond L + M = 15 we generally 
find an excellent agreement of the Pade approximants with 
the exact NCSM result for the same model space. This can be 
understood in terms of the Pade conjecture, which postulates 
the existence of a convergent subsequence of diagonal Pade 
approximants Nonetheless, there are individual 

approximants that show larger deviations from the exact re- 
sult, e.g. for d — 3, 4, 5 or 8. However, we observe deviations 
only for non-diagonal approximants, which are not covered by 
the Pade conjecture. In principle outliers are possible also for 
diagonal approximants, since only a subsequence of approxi- 
mants is expected to converge. 

The efficiency of the Pade resummation in recovering a ro- 
bust and accurate approximation from the divergent perturba- 
tion series is impressive, in particular for d = 3 and 4, which 
are the extreme cases of monotonous divergence. We observe 
that Pade approximants with at least L + Mx 10 are needed to 
obtain a quantitative agreement with the NCSM result. In turn 
this shows that the information contained in the high-order en- 
ergy corrections, which are responsible for the break-down of 
the power series, is indispensable to obtain stable and accurate 
results from the sequence of Pade approximants. 

In Fig. [3] we show the excitation spectrum for the positive 
parity states of 6 Li again for SRG parameters a — 0.04 fm 4 
(a) and a — 0.16 fm 4 (b). The columns from left to right 
represent the experimental spectrum, the exact NCSM results, 
the results from the Pade approximation and the spectrum ob- 
tained by truncating the power series at second, third, 4th, and 
8th order. We extract an averaged Pade result from the data 
shown in Fig. [2] in the following way: As discussed earlier, 
the diagonal as well as the sub- and super-diagonal Pade ap- 
proximants [L/M] for L + M > 15 are very stable, except for 
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FIG. 2: (color online) Pade approximation for energy of the ten states corresponding to the n = HO subspace of 6 Li. Shown are the 
diagonal [L/L] (•), super-diagonal [L/L + 1] (♦) and sub-diagonal [L/L - 1] (a) Pade approximants. The dashed horizontal lines represent 
the exact energies obtained by a NCSM calculation. We used a yV max = 8 model space with HO frequency 20 MeV and the SRG parameter is 
0.04fm 4 (A * 2.24 fnT 1 ). 
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FIG. 3: (color online) Low-lying spectrum of 6 Li computed in an 
A'max = 8 model space for hCl = 20 MeV, using SRG transformed 
N 3 LO two-body interactions with flow parameter 0.04fm 4 (A « 
2.24fnT 1 ) in (a) and 0.16fm 4 (A w 1.58 fm" 1 ) in (b). Shown are 
from the left the experimental observed energies, the NCSM results, 
the Pade resummed results (see text for details) and results from the 
truncated power series of DMBPT at order p = 2, 3, 4 and 8. Exper- 
imental values taken from Ref. |22T| . 



extremely few outliers. Therefore, we compute the average 
of all approximants for L + M > 15 excluding those approxi- 
mants which deviate by more than 0.5 MeV from the average 
of the remaining set. We use the standard deviation of the ap- 
proximants as a measure for the uncertainty of this average 
of the approximants. In Fig.[3]the bars representing the Pade 
results show an additional orange band representing this un- 
certainty. Overall, the uncertainty bands are very small and 
for the first 5 excited states we find remarkably good agree- 
ment between the Pade resummed excitation energy and the 



exact results obtained in the NCSM. For the four high-lying 
states we observe a small deviation of the Pade resummed re- 
sult from the NCSM values, reaching about 500 keV for the 
highest state. 

In contrast, the low-order DMBPT results do not provide 
a stable and reliable approximation for the exact eigenvalues. 
The excitation energies from second-order DMBPT are gener- 
ally too large and the deviations from the exact NCSM results 
increase with increasing excitation energy. The third-order 
contribution typically lowers the excitation energy. In a few 
cases the exact eigenvalue appears between the second- and 
third-order estimate, however, in other cases both low-order 
estimates are still above the exact NCSM eigenvalue. The in- 
clusion of the fourth-order contribution typically improves the 
results, but going to still higher orders destroys the agreement 
again. All these statements hold also for the soft potential with 
o- = 0.16fm 4 , see Fig. Hb). 

A similar picture emerges for the negative-parity spectrum 
of 7 Li depicted in Fig. [4] again using an N m - dK = 8 model 
space for the NCSM and the DMBPT calculations. Again, 
we find excellent agreement of the Pade-resummed energies 
with the exact NCSM spectrum for all states, with very stable 
Pade approximations giving rise to very small uncertainties. 
The only exception is the highest excited state computed with 
the a - 0.16 fm 4 interaction, that shows about 500 keV de- 
viation compared to the NCSM result. In contrast, the low- 
order DMBPT results show sizable deviations from the exact 
NCSM energies and change substantially from order to or- 
der. Again the second-order result typically overestimates the 
excitation energy and the third-order contribution lowers the 
excitation energy and often underestimates the excitation en- 
ergy. With increasing order the changes become less coherent 
and for high orders, beyond the order p — 8 the onset of diver- 
gence of the DMBPT series destroys the excitation spectrum 
completely (not shown in Figs. [3] and HJ. 

Finally, we note that though the general structure of the 
spectrum obtained in NCSM or Pade-resummed DMBPT is 
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FIG. 4: (color online) Excitation energies of the nine energetically 
lowest states of the 7 Li spectrum computed in an N mm = 8 model 
space for HQ = 20MeV, using SRG transformed N 3 LO two-body 
interactions with flow parameter 0.04 fm 4 (A » 2.24 far 1 ) in (a) and 
0.16fm 4 (A « 1.58 fm -1 ) in (b). Shown are from the left the ex- 
perimental observed energies, the NCSM results, the Pade resumed 
results (see text for details) and results from the truncated power se- 
ries of DMBPT at order p = 2, 3, 4 and 8. Experimental values taken 
from Ref. 121 . 

in agreement with experiment for the low-lying states, the ex- 
citation energies are systematically to high. This hints at defi- 
ciencies of the SRG-evolved two-body interaction used here, 
e.g. the lack of a three-nucleon interaction. The inclusion 
of three-body forces is straight-forward, because the recursive 
formulas of section [El] remain unchanged since the formalism 
is developed in terms of A-body Slater determinants. How- 
ever, this is beyond the scope of this paper. 



the way for first direct applications of DMBPT to open-shell 
nuclei and excited states, where one has to deal with degen- 
eracies by construction. 

We demonstrate the application of this formalism for the 
ground and excited states of 6 Li and 7 Li using an SRG- 
transformed NN interaction from chiral EFT at N 3 LO. In or- 
der to allow for a direct comparison with NCSM calculations 
for the same model space, we use the harmonic oscillator as 
unperturbed Hamiltonian. We find that the perturbation series 
itself is divergent for all states considered although the inter- 
actions used here are extremely soft and are sometimes termed 
'perturbative'. This divergence of order-by-order MBPT was 
observed for ground states of closed-shell systems already in 
Ref. 0]. To overcome this problem we resum the perturba- 
tion series using Pade approximants, i.e., we use the energy 
corrections from DMBPT to construct a rational function in- 
stead of a simple power series to approximate the energy. We 
find that Pade approximants using only low-order DMBPT re- 
sults yield no improvement compared to low-order DMBPT 
results. However, if we include DMBPT information from 
high orders, i.e., 15th to 30th order, the diagonal and neighbor- 
diagonal Pade approximants are in excellent agreement with 
exact NCSM calculations for the same model space. In con- 
trast, low-order DMBPT results are clearly not sufficient to 
provide a quantitative and reliable approximation to the exact 
NCSM excitation energies. 

A number of interesting topics remain, which we will ad- 
dress in the future: We can use the formalism of degenerate 
many-body perturbation theory to study binding-energy sys- 
tematics also for heavier open-shell nuclei, eventually also in- 
cluding three-body forces. Moreover, the convergence pattern 
of the DMBPT series will be affected by the partitioning of 
the Hamiltonians, i.e. by the choice of the unperturbed basis. 
There are hints from preliminary studies in the nondegenerate 
case that, e.g., a Hartree-Fock basis can improve the conver- 
gence behavior of the corresponding power series. 



V. CONCLUSIONS 

We have discussed degenerate Rayleigh-Schrodinger 
many-body perturbation theory for the description of ground 
and excited states of light nuclei with realistic Hamiltonians. 
We have derived an efficient, recursive formulation for the en- 
ergy and state corrections in DMBPT that allows us to evalu- 
ate the perturbation series order-by-order up to very high or- 
ders, typically up to p = 30. These formal developments pave 
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